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We study the nature of the zero-temperature phase transition between a d-wave superconductor 
and a Mott insulator in two dimensions. In this "quantum confinement transition" , spin and charge 
are confined to form the electron in the Mott insulator. Within a dual formulation, direct transi- 
tions from d-wave superconductors at half-filling to insulators with spin-Peierls (as well as other) 
order emerge naturally. The possibility of striped superconductors is also discussed within the dual 
formulation. The transition is described by nodal fermions and bosonic vortices, interacting via 
a long-ranged statistical interaction modeled by two, coupled Chern-Simons gauge fields, and the 
critical properties of this model are discussed. 



I. INTRODUCTION 

In recent years, due to remarkable experimental 
progress the cuprate superconductors have revealed 
a host of mysterious phases as their chemical doping is 
varied. Indeed, it would seem as though these materi- 
als exhibit many of the wide range of behaviors possible 
for low-dimensional, highly-correlated electron systems. 
Centrally located within the phase diagram and adja- 
cent to many of these puzzling regions is the d-wave su- 
perconductor. Beginning in this well-understood phase, 
one may develop theoretical descriptions of other, non- 
superconducting phases. Of particular interest are the 
T = quantum phases, both in the very underdoped 
and heavily overdoped regimes. The schematic situation 
is shown in Fig.|l|. 
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FIG. 1. Schematic phase diagram for the high Tc cuprates. 

When describing a 2-dimensional superconductor, 
topological defects in the Cooper pair wavefunction (BCS 
vortices) are of particular importance. Being bosonic, 
once they proliferate, they condense at T = 0, destroy- 
ing superconductivity 0] . In this way, a description 



of quantum phases with strong pairing correlations but 
lacking the phase coherence that is superconductivity 
emerge quite naturally as vortex condensates . If the 
superconductor is d-wave, there is the additional compli- 
cation of low-energy quasiparticles. As recently empha- 
sized [Q, there is a statistical interaction between these 
spin-carrying quasiparticles and the vortices, making the 
resulting theory strongly-interacting. 

Singlet-paired superconductors can be recast in a spin- 
charge separated form : the condensate carries all the 
charge but no spin, while the quasiparticles are electri- 
cally neutral with spin 1/2. Most other well-understood 
phases of electrons (such as the Fermi liquid) are spin- 
charge confined. It was recently argued that many 
puzzling aspects of the cuprate phase diagram could be 
understood in terms of the fractionalization and confine- 
ment of electrons. In this approach, the regions con- 
taining the pseudo-gap (and superconducting) phase are 
characterized by the presence of spin-charge separation 
(electron fractionalization) while the heavily over- 
doped regions are spin-charge confined. Between the two, 
a quantum confinement transition might cause critical 
behavior over wide regions of the high Tc phase diagram. 
A key feature of the cuprates is the close proximity be- 
tween a d-wave superconducting phase and a Mott in- 
sulating phase. Here, we work at half-filling and look 
at direct transitions between these two phases, viewing 
them as confinement transitions. We seek to answer two 
broad questions regarding the nature of such a transition. 

First, in terms of phenomenology, what sorts of states 
might we find when separate spin- and charge- carrying 
excitations are confined to form electrons? As we shall 
see, a remarkable feature of superconductivity with one 
electron per unit cell is that in the dual theory, the vor- 
tices are fully frustrated. When the vortices proliferate 
and condense, this frustration leads to the existence of 
multiple vortex condensates which break spatial symme- 
tries. In particular, we find vortex condensates which 
destroy superconductivity (and at half-filling, describe 
Mott insulators) as well as vortex condensates which pre- 
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serve superconductivity. As a consequence of the vortex 
frustration, we find direct transitions from superconduct- 
ing states to insulating states which spontaneously break 
rotational and/or translational symmetries, as well as 
the existence of superconducting states which have non- 
trivial spatial structure. Although we work at exactly 
one electron per unit cell, where the vortex theory is 
fully frustrated, it is hoped that even away from half- 
filling the qualitative features of our results will remain 
valid, in particular, the tendency toward spatial modu- 
lation near half-filling. In general, we hope that our ex- 
plorations of frustrated vortex systems can yield insights 
into quantum phases of electrons which are complicated 
by the presence of competing interactions. 

Second, as a specific example, we look at the critical 
properties of the the confinement transition between a 
spatially-modulated d-wave superconductor and a Mott 
insulator with the same broken translational symmetry. 
Characterized by the presence of long-ranged statistical 
interactions, which affect the confinement of spin and 
charge, this quantum critical point should have inter- 
esting universal properties. Within a special region of 
parameter space, we explore this transition analytically 
using renormalization group (RG) methods. 

Before we begin to address these questions, we first lay 
out the basics of the model under consideration. This 
model was introduced in Ref. [Q and many of its justi- 
fications and consequences can be found therein. Here, 
we provide only a whirlwind tour of its derivation and 
usefulness. 

II. THE MODEL: Z2 GAUGE THEORY 

We begin by formally writing the electron creation op- 
erator as a product of two operators, one of which carries 
the spin of the electron and the other, the charge. These 
operators are defined with the singlet-paired supercon- 
ductor in mind. If we write the Cooper pair creation 
operator as e*'^'', we construct our spinless charge e bo- 
son (called a "chargon") from the Cooper pair as: 

6t =s,e'^'-/' = e'*-,s, = ±1. (1) 

The chargon is "half a Cooper pair" in the sense that the 
square of bj. creates a Cooper pair. The neutral spin one- 
half particle (called a "spinon" ) is obtained by removing 
the charge from the electron: 

fra ~ ^rcli^- (2) 

As we shall see shortly, this spinon can be thought of as 
a neutralized BCS quasiparticle. With these definitions, 
we may perform a change of variables in a suitable Hamil- 
tonian describing electrons and Cooper pairs, resulting in 
a theory of chargons and spinous. However, the Hilbert 
space of chargons and spinons is much larger than that of 



electrons; for instance, the state with a single spinon but 
no chargons can be written down, but this state is un- 
physical and should be removed from the working Hilbert 
space. In other words, we may make this change of vari- 
ables only if we additionally impose a constraint that the 
sum of the number of chargons, Nr (canonically conju- 
gate to the chargon phase, [(j>r,Nr'] = iSrr' ), and the 
number of spinons, pr — fra fra, on each site is an even 
integer: 

This constraint can be implemented within a Euclidean 
path integral representation, resulting in a theory of 
spinons and chargons coupled to a Z2 gauge field [||. 
It should be noted that the constraint used here is not 
the same as Gutzwiller projection, and does not disallow 
doubly-occupied sites. 

For an odd number of electrons per unit cell and d- 
wave pairing correlations, the appropriate action in the 
Z2 gauge theory is: 

S^Sc + Ss + Sb, (4) 
S,^-t,J2'^^^ib:bj+h.c.), (5) 

<ij> 

-E7./^. (6) 

i 

where i and j label sites on a cubic space-time lattice. 
The Ising gauge field minimally coupled to the chargons 
and spinons, aij, can take values ±1, and Sb is a Berry's 
phase term. 

One may arrive at this action by making the above- 
mentioned change of variables in a Hubbard-type Hamil- 
tonian, as described in Ref. . Alternatively, this model 
can be taken as a starting point for describing systems 
with local singlet pairing correlations as well as Mott in- 
sulating tendencies. To exhibit the reasonableness of this 
model, consider the limits of infinite and vanishing tc- For 
tc — > 00, the bosonic chargons will condense and the Z2 
gauge field will become frozen with cr.y = 1, which frees 
the spinons. This phase is simply the d-wave supercon- 
ductor. The action reduces to S = Ss, which is just the 
Bogoliubov-deCennes action, with the spinons becoming 
the BCS d-wave quasiparticles. In the opposite limit, 
tc — > 0, the chargons are gapped into an insulating state. 
At tc — 0, the chargons may be trivially integrated out. 
The remaining action is just S = Ss + Sb- It is shown 
in Ref. [Q that the partition function for this remaining 
spin theory is formally equivalent to that of the Heisen- 
berg antiferromagnetic spin model. Therefore, we see 
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the attractiveness of this model for the cuprate system, 
which also exhibits both superconductivity and antifer- 
romagnetism. Many other additional properties of this 
action between these two limits are elucidated in Ref. , 
in particular, the presence of both spin-charge confined 
and deconfined phases |8|. 

The charge sector in Eqn. ^ is described in terms of 
the bosonic chargons, minimally coupled to a Z2 gauge 
field. In two spatial dimensions, vortices in the boson 
many-body wavefunction are point-like. This allows for 
a particularly elegant dual description where the vortex 
rather than the chargon is the fundamental degree of free- 
dom. In this duality, the condensate of chargons (the su- 
perconductor) is the vacuum of vortices; the condensate 
of vortices is an electronic insulator, where the chargons 
are gapped. Within the vortex theory, the superconduc- 
tor is trivial (being just the vacuum) and is therefore a 
good place to plant our feet. From this vantage, we look 
out of the superconductor at the neighboring insulating 
phases. The duality transformation, on the lattice, in the 
presence of the Z2 gauge field, has been explicitly imple- 
mented in Ref. [Q . The full resulting action at half-filling 
is: 

S ~ Ss + Sv + Sa + ScSi (8) 

= -t, J2 f^vCos{9, - dj + ^), (9) 

<ij> 

□ 

□ 

The spinon action Sg is unchanged. Here, e*^' creates 
an vortex and the flux of the U(l) gauge field, Uij, 
is the total electrical current. In particular, a flux of 27r 
through a spatial plaquette represents a charge of e. The 
terms Sy and Sa together form the usual dual vortex rep- 
resentation for charge 2e Cooper pairs except that here, 
the vortices are minimally coupled to the additional (Z2) 
gauge field /iy = ±1. The BCS vortex and the spinon 
are relative semions; upon circling a vortex, the spinon 
wavefunction picks up a minus sign. The term Scs is 
the Z2 analog of a Chern-Simons term for the two Z2 
gauge fields and mediates this statistical vortex-spinon 
interaction. The spinous "see" a Z2 flux J^g cr = —1 at- 
tached to each ^ vortex, while the vortices see a flux of 
Tin = (^1) ■ This flux attachment may be familiar to 
many in the context of the Quantum Hall Effect, where 
the gauge fields involved are for the U(l) group. Because 
of the anomalous "//" terms in the action, spinon num- 
ber is not conserved, and the usual Chern-Simons term 
cannot be used. 

In the superconducting state, we are in the vacuum of 
vortices. The spinous see no flux and are free to prop- 
agate independently of the chargons. However, when 



single vortices condense, the long-range statistical inter- 
action between the BCS vortex and the spinon drives 
spin-charge confinement. In the language of Ref. Q, the 
condensation of hc/2e vortices is accompanied by a con- 
densation of the visons (vortices in the Ising field, a), 
leading to a confined phase of electrons. 

We wish here to explore in some detail the nature of 
this confinement transition, where the freely-propagating 
spin and charge excitations are "glued together" to form 
the electron. Aspects and implications of this quantum 
critical point pertaining to the high Tc phase diagram 
have been introduced in Ref. |^. First, we will use Lan- 
dau theory to find phases related to the d-wave super- 
conductor by a second-order phase transition. Then, we 
will consider a special case where we recover a U(l) sym- 
metry for the spinous and will use quantum field theory 
methods to extract some analytic critical properties of 
the transition between deconfined and confined phases. 

III. DUAL VORTEX THEORY AT HALF-FILLING 

Concentrating on the vortices for the time being and 
working at half-filling, the dual theory for the charge sec- 
tor becomes 

£y = -ty cos{e^ - Oj + ^) + i |V X a - 27r'r|^ . (12) 

To obtain a low-energy effective theory, we work with a 
"soft-spin" model where the vortex creation operator e'^ 
is replaced by a complex field $. In the interest of ex- 
ploring the simplest case we set the charge per unit cell 
to be exactly e. In the dual theory, this corresponds to 
setting 

((V X a)f) = 27r. (13) 

In this section, we drop fluctuations of the gauge field a, 
and consider a Landau mean-field approach. This is jus- 
tified when the on-site repulsion between the electrons, 
U, is large. The vortices now see exactly (V x |)f = tt 
flux per spatial plaquette, and we are left with the two- 
dimensional fully-frustrated quantum XY model: 

J dr i ^ 2 - J2 ir>-'(*;-*f' + c.c.) 

+X][™'i$H'+"(i$rn']|, (14) 

where r labels sites on the 2d square lattice dual to the 
original electron lattice and the sign of t^^r' around a 
plaquette is —1. The sites of the dual lattice are at the 
centers of the plaquettes of the original lattice, and in 
units of the lattice constant (a — 1), r = {x,y) with x 
and y integers. 
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We proceed, following closely the work of others on 
the fully frustrated quantum Ising model by choos- 
ing the gauge (to be used in the remainder of this paper) 
seen in Figure |[ We may diagonalize the kinetic piece 
of this action to find two low-energy modes, living at 
{kx,ky) = (0,0) and (tt, 0), respectively. In real space, 
these (unnormalized) eigenvectors are: 



X°-=(l + V2) 



(l + \/2) 



(x,y integers). 



(1) 



(15) 
(16) 



(2) 



FIG. 2. Representation of the fully frustrated 2d lattice. 
Dark lines show the location of negative or "antiferromag- 
netic" bonds {try < 0). The two low-energy modes Xr and 
Xr are shown in (1) and (2), respectively. Long arrows have 
length 2 + ^/2 and short arrows have length y/2. Frustrated 
or "unhappy" bonds are marked with a slash. 



For the purpose of characterizing the low-energy be- 
havior of this vortex system, we consider fields which are 
linear combinations of these two low energy modes. 



$(f, t) «'o(f, r)x°- + «'^(f, t)x2 . 



(17) 



We now have two complex fields, 'i'o{f) and '^'^^{f), which 
describe the low-energy configurations of our vortex sys- 
tem. The phase transitions of the system can be ex- 
plored within Ginzburg-Landau theory. The Ginzburg- 
Landau Hamiltonian for the two- vortex system must pre- 
serve all the symmetries of the original lattice Hamilto- 
nian, namely: discrete x and y translations, rotations 
by ^, and the vortex U(l) symmetry ($ e*"<i>), as 
well as hermiticity. In terms of our two complex vor- 
tex fields, these symmetry transformations take a simpler 
form when expressed in terms of the fields 



as follows: 



Ty 



V2, 



(18) 



(19) 



Ri 



-101, 



C/(l) : cba~^ e'^cba (for a=l & 2). 
Allowed terms for the action include 



(20) 

(21) 
(22) 



(I) 
{II) 
(III) 



(|01P)" + (|02P)", 
(I01PI02P)", 
[(0t02)* + (010^)1" 



(with arbitrary positive integer, n), and combinations of 
these terms. Expanding in powers of the fields, we take 
as our Landau-Ginzburg action: 



SlG = f d^xdT J2 [\Ma 



2 + r|(i„,n 



+U^{llaWY + vMrYW - V^W^C^^r + /i.e.], (23) 

where r has been rescaled to set the vortex velocity 
— \. The terms labeled by M4 and W4 are the only 
allowed quartic terms, and are invariant under indepen- 
dent U(l) transformations on 0i and 02- We have kept 
the f8 term because it is the lowest-order term which 
breaks this symmetry down to the global J7(l) of Eqn. 
p^ . This model will be employed to construct a descrip- 
tion of various phases proximate to the d-wave supercon- 
ductor within mean field theory. 

We wish to characterize the various states of this vortex 
system. It is important to emphasize at this point that 
not all vortex condensates destroy superconductivity. Su- 
perconductivity is destroyed when the dual U(l) symme- 
try of the vortex theory (Eqn. p^ ) is broken. Therefore, 
it is possible to have non-trivial vortex condensates which 
are superconducting. This leads to two scenarios for the 
superconductor-insulator transition at half-filling. First, 
we may consider superconductors which arc described by 
a vacuum of vortices; superconductivity is then destroyed 
when single vortices proliferate and condense (in a way 
which breaks the dual U(l)). Alternatively, the super- 
conducting state could itself be a U(l)-preserving vor- 
tex condensate which then undergoes a transition which 
breaks the dual U(l), killing superconductivity. 

In the following sections, we explore the phases of our 
dual vortex model using the Landau-Ginzburg action of 
Eqn. 23. Due to the frustration of the vortex theory 
with one electron per unit cell, the vortex condensates 
will break lattice symmetries. Some of these spatially- 
ordered states are superconductors and some are insu- 
lators. We will begin by describing the possible super- 
conducting states within the dual theory (including, a 
striped superconductor), and then move on to a descrip- 
tion of the insulating states. Ignoring charge fiuctua- 
tions in the superconducting states (as we have in arriv- 
ing at Eqn. p3) is not justified, and a good description of 
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these states would require putting the charge fluctuations 
back in. Here, we content ourselves to characterizing the 
phases of our vortex system by their broken symmetries. 
We conclude with a summary of the possible transitions 
from superconductor to insulator within this mean field 
theory. 



A. Superconductors 

1. Vortex Vacuums 

The simplest superconducting phase is just the vortex 
vacuum. This is the standard BCS d-wave superconduc- 
tor. Destruction of superconductivity occurs when sin- 
gle vortices proliferate out of the vacuum and condense, 
breaking the dual U(l) symmetry. The effective action 
for this transition is Eqn. 

2. Paired Vortex Condensates 

Condensation of single hc/2e vortices necessarily 
breaks the dual U(l) symmetry (Eqn. |2^), killing super- 
conductivity. However, when pairs of vortices condense, 
the U(l) can be preserved. Consider the paired vortex 
condensate: 



7^0, 



0. 



(24) 



We see that in this condensate, the dual U(l) is pre- 
served, and the state is characterized by the phase of 
the condensate (setting the amplitude |(02'/'i)| = 1 for 
simplicity) , 



1)291/ 



,ix)-d2ix). 



(25) 
(26) 



Here, 9i and 62 are the phases of 0i and ^2, respec- 
tively, and are still free to fluctuate. Only the combi- 
nation 6 = 61 — 02 is uniform, reflecting the fact that 
the dual U(l) symmetry is preserved (i.e., (fii and (j>2 are 
uncondensed) . The only term in the Landau-Ginzburg 
action which depends on 9 is the vg term, giving: 



(27) 



Sy = —vs / d xdTCOs{4:9). 



We see that the ground state depends on the sign of «§: 



vs>0 



n 



2' 



„ TT TT 



(28) 
(29) 



with n an integer. 

The spatial symmetries in EqnsjTg|-pT|, written in terms 
of the relative phase 9, are: 



Ty 



(30) 
(31) 
(32) 



From this we can see that the vortex condensate fa- 
vored by ug > breaks the lattice rotational symme- 
try and one of the two translational symmetries. We 
therefore associate this condensate with a stripe-type or- 
dering: a striped superconductor. This state is particu- 
larly interesting given recent experimental results which 
suggest possible stripes in the superconducting state of 
La2-xSrj;Cu04 [|l0|. The ground state for vg < breaks 
all of the lattice symmetries; we identify this state with 
a "plaquette" order which will be made more explicit in 
upcoming sections when we discuss the insulating states 
of the vortex system. For now, we emphasize the possi- 
bility of spatially-ordered superconducting states which 
emerge quite naturally within our dual vortex descrip- 
tion. 

Still working in the dual description, these striped and 
plaquette superconductors are described by an effective 
theory of one vortex species, since the paired condensa- 
tion has locked the two original vortices together: the 
vortex phases 9i{x) and 92{x) — 9i{x) — 9 still fluctu- 
ate within the superconducting phase, but not indepen- 
dently. When the remaining phase 9i becomes constant 
over the sample, the dual U(l) is broken, and supercon- 
ductivity is destroyed. Therefore, for these spatially or- 
dered superconductors, Slg (Eqn. |2^) reduces to: 



Sy 



d^xdr [|a^</.ip-f r|0i|2 + u(|0i|2)21 



(33) 



It is worth noting that we have gone from a theory of a 
single fully-frustrated vortex to a theory of a single un- 
frustrated vortex via a theory of two vortices. This is 
possible because in a striped or plaquette superconduc- 
tor, the unit cell is doubled. If one started from scratch 
in constructing a dual theory of these striped (plaquette) 
superconductors, the vortices would see 2tt rather than tt 
flux per (doubled) unit cell and there would be only one 
low-energy mode. 



B. Confined Insulators 



When single vortices condense at half-filling, we move 
from the d-wave superconductor into a confined insula- 
tor. Within our dual formulation, these insulators are 
described by condensates which break the dual U(l) sym- 
metry of Eqn. In the case of superconductors which 
are vortex vacuums, because we have two vortex species, 
there are many ways to do this and therefore many pos- 
sible single vortex condensates. We will see that these 
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different vortex condensates correspond to different insu- 
lating states of electrons. We return to the case of the 
striped and plaquette superconductors after first enumer- 
ating the insulating states at the mean field level, using 
the action of Eqn. |2^ . 

The most general U(l)-breaking vortex condensate is: 



(34) 
(35) 



where |(02)|, Oi, and 62 are all fixed real numbers. 

Within our dual Landau-Ginzburg model, condensing the 
vortices corresponds to setting r < and U4 > 0. The 
signs of U4 and then determine the ground state. For 
W4 < 0, both vortex species acquire a non-zero amplitude 
\{'Pi) \ = 1(02)1 7^ and their relative phase 6*12 =61 — 62 
is determined by the sign of vg. On the other hand, if 
V4 > 0, the ground states are condensates of either i/ii or 
<j)2 and the sign of us is irrelevant. Each of these conden- 
sates will correspond to a different insulating state of the 
electron system. We consider each case in turn. 



1. K0i)l = K<^2)l/o 



of the system in this case are characterized by "stripes" 
of energy-density as shown in Figure 0b . 



(a) 









\ 
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\ 


\ 









(b) 



These condensates are favored by U4 < 0, and the rel- 
ative phase {9 1 — 62) is determined by the sign of vg. 
Taking the magnitudes \{(f>i)\ = \{4'2)\ = 1, this term in 
the action can be rewritten as: 



V8[{(l)*i(l)2)^ + h.c] = -vs cos(46li2). 



(36) 



In terms of this relative phase, the spatial symmetries are 
given by Eqns. pQ-p2 with the replacement 6 ^ 612- 



FIG. 3. (a)The vortex condensates with 612 = f , which 
correspond to $(r-) = + Xr)- Frustrated bonds are 

slashed, (b) Schematic of the four "striped" states corre- 
sponding to = with higher energy (frustrated) bonds 
slashed. 



This class of condensates is preferred hy vg > 0. There 
are four general states, corresponding to each of the possi- 
ble values of n. We see by the symmetry transformations 
in Eqns. ^0[]3^ that each of these states breaks the lattice 
rotational symmetry as well as breaking one of the two 
translational symmetries while leaving the other intact. 
On these grounds alone, we could guess that these states 
correspond to "stripe-like" phases. To be more concrete, 
we may go back to our real-space representation for the 
vortex field <i>(r) in EqnjlJ and draw real-space pictures 
of these lattice states. The values of the fields at various 
points will be gauge dependent, but the location of frus- 
trated bonds (which are places of higher energy density) 
is gauge- independent and therefore a good way to charac- 
terize the state of the system. This is shown for the case 
di2 = J, as an example, in Figure Investigations of 
this sort lead us to conclude that the four ground states 



We now turn to a characterization of this system in 
terms of the electron degrees of freedom. Because we 
have broken the dual U(l) symmetry of the vortices and 
we are at half-filling, these states will be Mott insula- 
tors. The charge degrees of freedom live on the plaque- 
ttes of the dual lattice and are fixed at one charge of e 
per dual plaquette. It has been suggested fll],^ that the 
frustrated bonds of the dual lattice should correspond to 
singlet bonds of the electron system, since one expects 
regions of higher energy along the links where the elec- 
trons spend most of their time. This relationship between 
the frustrated bonds on the dual lattice and the singlet 
bonds on the original lattice is illustrated in Figure |^. 
These "striped" vortex phases then correspond to spin- 
Peierls (or "bond density wave") order in the insulating 
electron system. 



6 



FIG. 4. Schematic of the relationship between frustrated 
bonds of the dual lattice (slashed) and links of the original 
lattice where the corresponding "singlet bonds" live (dashed) . 



FIG. 6. Schematic of states corresponding to the vs < 
vortex condensates. Relatively unhappy bonds of the dual 
lattice are slashed, and links of the original lattice where the 
singlet-type bonds live are dashed. We associate this struc- 
ture with a "plaquette RVB" state of the electrons. 



i) / 0, (02) = or {^^) = 0, (^a) / 



These condensates are favored by < 0. Here, how- 
ever, each ground state breaks all of the discrete lattice 
symmetries. We proceed as above and obtain character- 
izations of these states in terms of the location of frus- 
trated bonds. The result is a plaquette-like structure, as 
seen in Figure ^ 



■^^ — ■ — ^ 



4-1 1 1 1 1 1 1 1 i :t 1 1 1 1 1 1 1 1 



I I I I I I Ut I I I I I 



These condensates are preferred in the case 'i;4 > 0. 
We may proceed as above in drawing real-space diagrams 
corresponding to these states. We find, as shown in Fig- 
ure ^ that these states have vortex currents around each 
plaquette, of alternating sign. 



^ — s — 



J 































FIG. 7. The vortex condensate with <3?r = Xr + *Xr • Since 
the field is complex, the magnitude is given by the length and 
the argument by the direction of the arrow at each site. The 
direction of vortex current is indicated on the bonds. 



FIG. 5. The vortex condensate with 6112 — j, which cor- 
responds to $r = 0.9xr + 0.4Xr • Locations where the field 
is zero are denoted by a dot. Relatively unhappy bonds are 
slashed. 



In order to interpret this state, we will have to put back 
in the spinons which have been ignored in the previous 
discussion. The vortex-spinon action is: 



In terms of the electron degrees of freedom, we would 
like to again interpret the frustrated bonds of the dual 
lattice as regions where singlet-type bonds of the elec- 
tron system reside. The plaquette-like structure of these 
vortex states may then correspond to a "plaquette RVB" 
state of the electron system, as shown in Figure ^. 



S — Sv + Ss + Scs, 
Sv = -ty y /i,j/,v cos( 



<i'j'> 



with (V X a) = 27r, 

5. = - V o,m,h!,^^!^^J,^- 



<ij> 



E 



fi fi 1 



(37) 
(38) 

(39) 
(40) 
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^c^^=E*i(i-n^)(i-'^^^)' (41) 

where i,j label sites on the original lattice and la- 
bel sites on the dual lattice. Looking at Sv, we see that 
the alternating vortex currents would like to induce com- 
pensating fluctuations in either the a^'j' or fii'ji fields. 
Allowing fluctuations of the gauge field a (which de- 
scribes charge fluctuations), and ignoring the coupling 
to the spinous, the alternating vortex currents would in- 
duce CDW order at wavevector {tt,7t). However, with a 
large on-site U, this state will be greatly suppressed. If 
we forbid charge fluctuations, we see that the alternating 
vortex currents will instead drive a mean fleld in the Z2 
gauge field: 

Hf,,,,, ^ {-ir ^ -1, (42) 
□' 

(where is the number of spinous in the dual plaquette 
denoted by □), which corresponds to one spinon per unit 
cell. Unlike the previously-considered vortex condensates 
(with V4 < 0), at the level of vortex mean field theory, 
this state has no broken translational symmetries. (How- 
ever, we cannot rule out the breaking of symmetries by 
the charge and spin fluctuations.) We note that a pos- 
sible candidate for this state which has one electron per 
until cell and uniform energy density is the antiferromag- 
net. 

C. Summary of Vortex Theory 

We have seen that our dual vortex theory describes 
both standard BCS and striped or plaquette d-wave su- 
perconductors as well as a host of confined insulating 
states. We now summarize the results of Landau theory 
for the transitions from the d-wave superconductor to the 
confined insulator at half-filling. 

We consider first the transition from a vacuum of vor- 
tices (a superconductor) to a U(l)-breaking condensate 
of vortices (an insulator). Within mean field theory, the 
nature of the insulating state is determined by the signs 
of the coupling constants in Eqn and we have the 
following possible direct transitions out of the symmetric 
d-wave superconductor: 

W4 < 0, ^8 > 0; dSC spin-Peierls, (43) 
V4 <0, vs < 0; dSC plaquette RVB, (44) 
W4 > 0; dSC uniform state of electrons. (45) 

One might hope to ascertain which of these insulating 
states is preferred close to a d-wave superconductor, in- 
cluding fluctuations beyond the mean field level, by con- 
sidering the fixed points of the action in Eqn.^. In par- 
ticular, we see that the sign of W4 determines whether we 
enter one of the states of broken translational symmetry 



(spin-Peierls or plaquette RVB), or the state with uni- 
form energy-density (possibly the antiferromagnet) . The 
work of Blagocva ||l^ on the theory of two-component 
complex fields with these (and other) couplings gives a 
stable fixed point at 114 < 0, to order (e = 4 — D, 
D ^ d + 1, in d spatial dimensions). This suggests that 
the transition dSC —f spin-Peierls would be preferred over 
dSC — *■ uniform state. This is tantalizing given the ex- 
perimental evidence for intervening "stripey" phases be- 
tween the superconducting and antifcrromagnetic phases 
in the cuprates |l3| . 

In the case of the striped and plaquette superconduc- 
tors, when the single vortex species in Eqn. |3^ condenses, 
superconductivity in these states is destroyed and we 
enter a confined insulating state. Because the relative 
phase 01 — 62 is already fixed within these superconduc- 
tors, we see from our above analysis of the insulating 
phases that the insulating state is pre-determined. The 
striped superconductor (with 9 = mr/2) enters the spin- 
Peierls insulator, and the plaquette superconductor (with 
6 = 7r/4 + n7r/2) enters the plaquette-RVB insulator. 
In other words, these spatially-ordered superconductors 
make transitions into insulating states with the same bro- 
ken spatial symmetries: 

vg > 0; striped SC spin-Peierls, (46) 
vg < 0; plaquette SC plaquette-RVB. (47) 

In the preceding section, we have considered states of 
electron systems at half-filling near a d-wave supercon- 
ductor within a dual formulation in terms of vortices. 
Each phase is characterized by a dual (vortex) order pa- 
rameter. At one electron per unit cell, the frustration of 
the vortex theory manifests itself in spontaneously bro- 
ken spatial symmetries. Exploiting the fact that the 
vortex order parameters break spatial symmetries has 
helped us identify these vortex phases with more famil- 
iar phases of electrons (such as the spin-Peierls state), as 
well as phases like the striped superconductor. The dual 
formulation shows us the enhanced chance for striped 
superconductors near half-filling. In the next section, we 
will add back in the spinous (and along with them, their 
long-ranged statistical interaction with the vortices), and 
extract information about the critical properties of the 
confinement transition using field theory methods. 

IV. CONFINEMENT TRANSITION 

Having explored the vortex sector of the theory with 
one electron per unit cell, we now wish to put the spinous 
back in and address the critical properties of the confine- 
ment transition. Because we will continue to work at 
half-filling, the confined states of electrons will be Mott 
insulators. While the theory of vortices and spinous cou- 
pled to Z2 gauge fields may in principle be numerically 
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accessible, the action suffers from the notorious fermion 
sign problem. Here, we discuss a special case which will 
turn out to be accessible to perturbative RG calculations. 

Focusing on the spinon Hamiltonian (and dropping the 
Zi gauge field for the time being): 

= - E + C.'(/.T/r'i + c.c.)], (48) 

<rr'> 

we choose the special case: 



with nearest-neighbor d-wave pairing amplitude: 



i'^ — ~f 
^r,r±y ~ 

Following Affleck jlj], we introduce the fields 



(49) 



(50) 



for y even. 



{~icry)e ^s'^y [ ) , for y odd, 



(51) 



(where ay is the usual Pauli matrix), the spinon Hamil- 
tonian becomes: 



i?s = - E Ur'idladr'a + h.c), 



with 



—t for y and y' even, 
t else. 



(52) 



(53) 



This is the Hamiltonian of fermions hopping in 2d in the 
presence of tt flux per plaquette. We have succeeded in 
finding a Hamiltonian for the spin sector which has a 
conserved fermion number. The original theory (Eqn.^) 
can now be written in terms of these d fermion fields, the 
chargons and the Z2 gauge field. Following a transfor- 
mation which can get rid of the Berry's phase term fl^ ] 
this theory can be modeled numerically with no fermion 
sign problem. Here, we instead proceed to a low-energy 
continuum Hamiltonian for the spin sector. To this end, 
we diagonalize Hs to find two Dirac points. These are the 
usual d-wave quasiparticle nodes at {k^, ky) — (±^, ±-|) 
except that due to the tt flux per plaquette, we have dou- 
bled the unit cell and halved the Brillouin zone; it now 
contains only two of these nodes, which we denote Ki and 
K2 ■ In terms of long- wavelength fields living at these two 
nodes. 



dja{x) ~ iJijaix^e^^^'"" + ^2ja{x)e 



(54) 



(where j = 1,2 labels the sublattice), the continuum 
Hamiltonian is: 



Hs = j cPx Vs?/>|^[Ti(-i9^) 4- T2{~'idy)]i:ia 

+Vsl/'L['r2(-»5a;) + Ti{-idy)]ll)2a, (55) 



where, 



r.^i=(r.-r.). 



(56) 
(57) 



Here, r acts in the sublattice space, and we have ro- 
tated the X and y axes at each node by 45°. Just as the 
Hamiltonian for the d fermions was diagonal in the spin 
label, so is this one, and we are left with a theory of four 
species of Dirac fermions. Note that the spinon charac- 
teristic velocity, Vs is isotropic in space because of our 
choice ts = \t\\. 

Defining Dirac matrices in 2-1-1 dimensions: 



at node Ki : at node K2 : 



70 

71 = T2, 

72 = -n, 



70 

71 = -Tl 
72 = T2, 



(58) 



(7m)^ = 7p ; {7m'7i^} = 2(5^^, (at each node), (59) 
we proceed to the Euclidean Lagrangian density: 



V' = "0^0 • 



(60) 
(61) 



The fields if^a have two components (corresponding to 
the sublattice label), and summation conventions on the 
number of species a S [1,4] (one for each spin at each of 
the two nodes) and the spatial dimension i e [1,2] are in 
use. 

We have succeeded in writing a low-energy effective 
theory for the spin sector which is just that of four species 
of two-component Dirac fermions. We may now write 
down a full low energy effective theory where, since the 
spinon and vortex sectors each display U(l) symmetries: 



(62) 
(63) 



we may implement the statistical spinon-vortex interac- 
tion using U(l) (rather than Z2) Chern-Simons gauge 
fields, A*^ and A'''. We proceed to a field theory model- 
ing the confinement transition between a spin-charge sep- 
arated (d-wave) superconductor and a spin-charge con- 
fined Mott insulator. For simplicity, we consider the vor- 
tex theory with only one species (Eqn.^Sj), and thereby 
consider transitions out of the striped/plaquette super- 
conductor given in Eqns. ^ - The low-energy effec- 
tive action is: 
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5* = 



Cy+£ 



cs 



(64) 



r 

C-cs - 
r 

with 



i,^^ - igP)i^a + ^^aM - ^91:At)^a, (65) 



|(9^-zg^^P+m2|0|2 + i,o( 
ae [l,iV = 4] , K=—-l, 



|2\2 



(66) 
(67) 
(68) 
(69) 



where k is a measure of the velocity anisotropy between 
vortices and spinons and will be treated as a perturba- 
tion. We have added the term C^s in the interest of 
including all possible relevant interactions. The Chern- 
Simons term causes a vortex taken around a spinon to 
acquire a phase of: 



o«9 



and likewise for a spinon after encircling a vortex: 



ig i, A^-dl 



(70) 



(71) 



so that the full statistical interaction is achieved when: 
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2Txa\ a 



(72) 



(where a is the so-called "statistics angle" and is equal to 
1/2 since the vortex and the spinon are relative semions). 
The theory as written neglects charge fluctuations, which 
is not justified within the superconducting phase. The 
full vortex theory would include an additional minimal 
coupling to a gauge field a |l^ . As seen in the dual XY 
model, this coupling causes runaway flows, and is prob- 
ably best modelled numerically. At this point, we leave 
out the gauge field a and its attendant problems, but we 
will revisit this question shortly. 

When the vortex Lagrangian is taken through criti- 
cality (m} < 0), the statistical interaction, mediated by 
the gauge fields Aj^ and A^, will drive spin-charge con- 
finement. Here, we seek the effect of these statistics on 
critical properties of the system. In particular, we wish 
to calculate /? functions for the couplings uq, vq, k, and 
(7, as well as the anomalous dimensions of the vortex and 
spinon fields. 

We work hi D — d -\- 1 = 3 dimensions (indeed, 
our Chern-Simons flux attachment is not well-defined in 
higher dimensions), and define dimensionless couplings: 



A ^uqKd. 

VoKd=3, 



(73) 
(74) 
(75) 



where factors of Kd = [2^-^TT^/^r{D/2)] have been 
put in for later convenience. The bare propagators in the 
Landau gauge are: 



termions: CTn, = — — r^-, 
vortices: Ct = —rr, 



gauge fields: S'q 



{A^A"^) = (A'^A^) = 0. 



(76) 
(77) 

(78) 
(79) 



(The fermion propagator is diagonal in the label a, so we 
have suppressed this index.) 

For the /3-functions we find, to lowest non-vanishing 
order (1 loop): 



du 

m 

dv 
dl 

dl 



= u- lOu^ + (^ + Cn)v^ + 
3 

— —Auv -I- • • • , 
= 0, 



dn 

1! = ' 



(80) 
(81) 
(82) 
(83) 



We expect that at higher orders, g will enter into ^ and 
^ non-trivially, but that g itself should not renormalize 
at any order, following the argument given by Semenoff 
et al. ||l^ . The one- loop RG equations for u and v have a 
stable solution at w = 0, so that the theory decouples into 
separate spinon and vortex theories. At this order, since 
the spinon and vortex sectors decouple, we may ignore 
the Chern-Simons gauge fields (effectively taking g — Q) 
and include the effects of charge fluctuations by using the 
full dual XY model for the vortex sector: 



r 



1(9^ - ieoa^ 



4' 

|2\2 



IV X al' 



(84) 



Recently, much work has gone into tackling the critical 
properties of the e ^ model , and we may use these 
results. 

To first order, then, we find a fixed line, parameterized 
by values of the statistics angle a (or equivalently, the 
coupling g). At lowest non- vanishing order, this line is 
given by: 



V* ~ 0, 

{gr = 9'-^, 



(85) 
(86) 
(87) 
(88) 
(89) 



where, by it^„^; and eduai, we mean the values of these 
couplings at the fixed point of the dual XY model. 

In order to see whether spinon-vortex velocity 
anisotropy grows, we need to take the j3 function for k to 
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its lowest non-vanishing order, which is two loops. The 
result is: 



dn 
H 



240 7r2 



(90) 



Since the system flows toward k = 0, it is legitimate to 
treat this term as a perturbation, and the theory becomes 
"relativistic" at the critical point. 

We proceed by calculating the anomalous dimensions 
of the spinon and vortex fields, to lowest order, near the 
critical point. To that end, we consider the self-energies: 



(91) 
(92) 
(93) 



Near the critical point, the anomalous dimensions are 
given by: 



G"^(fc) 
G^(fc) 



1 



|/c|2-'J0^ 



|A:|2-w 



(94) 
(95) 



(up to additive constants). Working at the fixed point 
{u,v,g'^) = {u* = u*^uav'"* = = ^^d calculat- 

ing the spinon and vortex self-energies to two loops in 3 
dimensions, we find: 



»?0 
'7^ 



4 (.9^) 



4n* 



3 IGtt^ 



■N - 



4\* 



3 167r2 



(96) 
(97) 



Since we are in the case with one vortex species, we 
may take the numerical results of Hove and Sudb0 p^ ] 
for the anomalous dimension of the vortex field in the full 
dual XY model in D = i: rj^ —0.24. After plugging in 

= 4, and = tt in our result, we find: 



r;^ ~ -0.24- i~ -0.57, 



48 



-0.02. 



(98) 
(99) 



These critical exponents may reveal themselves in 
many quantities. In particular, the spectral function as 
probed by ARPES and the spin-spin correlations probed 
by NMR or neutron scattering. Within our theory, the 
low-energy electron correlator decouples into chargon and 
spinon pieces for 0: 



z{x)cHO)) = {h{x)hHQ)){f{x)f\0)). 



(100) 



These correlators will exhibit anomalous dimensions T/f, 
and "qf, which can be expanded perturbatively around 
5* =0: 



rib=VxY + C^;;\g*Y + 



rjf ^1]^^ Cf\g*y -I , 



(2) 

/ 



(101) 
(102) 

-0.03. The anoma- 



where we have calculated C 

lous dimension for the 3d XY model (appropriate for 
one vortex species) has been calculated by Hasenbusch 
and Torok using Monte Carlo methods ||l^; they find 
rjXY — 0.038. The anomalous dimension will also en- 
ter into the spin-spin correlation function. Within our 
model, it looks as though vertex correction diagrams will 
not contribute as much near the critical point as the di- 
rect [G^]^ term. 



V. CONCLUSIONS 

In this paper, we have used a gauge theory of strongly- 
interacting electrons to explore the regions near the su- 
perconducting state in the high- Tc cuprates. This gauge 
theory exhibits spin-charge separated and spin-charge 
confined phases. We have seen that the presence of one 
electron per unit cell has profound implications for the 
regions near the superconducting state. Within a dual 
description, half-filling of electrons corresponds to fully- 
frustrated vortices, leading to a spontaneous breaking of 
translational symmetries in the electron system. From 
this, we have seen the possibility of striped superconduc- 
tivity as well as a host of confined insulators descending 
from d-wave superconducting phases. We have then used 
Chern-Simons methods to calculate lowest-order critical 
properties of the confinement transition between these 
phases. Because we have worked at half-filling of elec- 
trons throughout, our results are of particular relevance 
to the undoped cuprate materials, which may be spin- 
charge confined. However, we also hope that the flavor 
of our results may be of interest in the heavily overdoped 
materials, where the conflnement of spin and charge may 
result in a Fermi liquid phase. 
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